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Abstract 

We establish the Plancherel-Rotach-type asymptotics around the largest zero 
(the soft edge asymptotics) for some classes of polynomials satisfying three-term 
recurrence relations with exponentially increasing coefficients. As special cases, 
our results include this type of asymptotics for g -1 -Hermite polynomials of Askey, 
Ismail and Masson, g-Laguerre polynomials, and the Stieltjes-Wigert polynomials. 
We also introduce a one parameter family of solutions to the g-difference equation 
of the Ramanujan function. 
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1 Introduction 

Using the explicit representations and g-identities (in particular, the g-binomial theo- 
rem), Ismail [7] derived the complete asymptotics expansion of g _1 -Hermite, g-Laguerre, 
and Stieltjes-Wigert polynomials near their respective largest zeros. One of the main 
features is that the Ramanujan function A q (z), called the g-Airy function in |7]|, appears 
very naturally in all the results. These asymptotics are called the soft edge asymptotics 
in random matrices. There are two other types of Plancherel-Rotach asymptotics, the 
bulk scaling asymptotics where you normalize by the largest zero but keep x in the 
oscillatory region, and the tail asymptotics which is asymptotics beyond the largest and 
smallest zeros [18J. Ismail and Zhang [IT] p2] derived these asymptotics and found 
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that the sine and cosine functions which appear in the bulk scaling asymptotics of the 
Hermite and Laguerre polynomial asymptotics [IB] are now replaced by theta functions. 

In this work, we will start with the more general three-term recurrence relations 
that include these classical g-orthogonal polynomials as special cases and derive the 
Plancherel-Rotach type asymptotics. A general characteristics of these recurrence re- 
lations is that the recurrence coefficients are growing exponentially with n, the degree 
of the polynomials. The orthogonality measures of these polynomials are not unique. 
Through careful scaling and transformation, these recurrence relations can be trans- 
ferred into more trackable ones on which we base our study of the convergence. It must 
be emphasized that we use the recurrence relation directly in contrast to the previous 
works, [7], [EE] [12], which use more detailed information, such as explicit formulas. 

The approach of deriving the asymptotic properties (e.g., ratio asymptotics, zero 
distributions, etc.) of orthogonal polynomials from assumptions on the asymptotic be- 
havior of recurrence coefficients can be traced back to the work of Poincare [16] and 
Blumenthal [3] . See the monograph in the Memoirs of AMS by Nevai [H] for a system- 
atic treatment of the class M(a, b) (the Nevai-Blumenthal class as it is referred to now, 
see, e.g., [13]) where the recurrence coefficients are assumed to be convergent. Starting 
with the assumptions on the recurrence coefficients, Van Assche [20] and Van Assche 
and Geronimo [22J obtained the Plancherel-Rotach asymptoics outside the oscillatory 
region (see also [21] and the references therein). Recently, Tulyakov |T9] proposed a 
new method for obtaining the Plancherel-Rotach type asymptotics when the recurrence 
coefficients are rational functions in n (indeed, he dealt with more general difference 
equations of order p). This method allowed him to obtain the global picture of the 
asymptotic behavior of the solutions, and as a demonstration, he applied the method to 
the Hermite and Meixner polynomials. 

In recent years the Riemann-Hilbert problem approach has been a powerful technique 
in determining asymptotics of orthogonal polynomials, see [4J for the description of the 
Riemann-Hilbert approach and [5] for the application to derive asymptotics of polynomi- 
als orthogonal with respect to exponential weights. So far the Riemann-Hilbert problem 
has not been successfully applied to g-orthogonal polynomials except in the case of the 
Stieltjes-Wigert polynomials when q = e _1 ^ 4n ) and n is the degree of the polynomial, 
see 0. 

Throughout this work, we will assume that < q < 1. 

Our model case is the g _1 -Hermite polynomials. Recall that the g _1 -Hermite polyno- 
mials {h n (x | q)} of Askey pQ and Ismail and Masson [10] satisfy the recurrence relation 

(1.1) 2xh n (x | q) = h n+1 (x | q) + q' n (l - q n )h n ^(x | q), 

with ho(x | q) — 1 and h\{x \ q) = 2x. For detailed properties of the g _1 -Hermite 
polynomials see [H Chapter 21]. Using a result in [9] it is easy to see that the largest 
zero of h n (x \ q) is less than (but close to) q~ n l 2 . Let 

(1-2) x n (t)= l -[q~ n /H-q n / 2 /t\. 
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Clearly, x n ±i(q ±1 / 2 t) = x n {t). Set x = x n (t) in fll.ip to get 

re/2" 



q 



-n/2 t 1 



, K(x n (t) | g) = h n+1 {x n+1 (q^H) | g) + g-"(l-g"')/i n _ 1 (x n _ 1 (g- 1 / 2 t) | q). 

With p„(t) = t~ n q nI / 2 h n (x n {t) \ q), the above recurrence becomes 

(1-3) (l - £\ Pn (t) = Pn + M /2 t) + l^fp^q-^H). 

If p n {t) — > f(t) uniformly on compact subsets of C \ {0}, as n — > oo, then f(t) will be 
analytic in C \ {0} and satisfy 

(1-4) f{t) = f(q 1/2 t) + j 2 f{q- l/ H) 

Since p n (t) is a polynomial in 1/t 2 , / will have the Laurent expansion: 

oo 

(1.5) f(t) = Y,U- 2n , With / = 1. 



n=0 



The substitution of / from (11.51) in (11 .4^ then equating coefficients of various powers of 
t implies that 

f{t)=A q {t-% 

where A q (z) is the Ramanujan function [17] which plays the role of Airy function, as 
pointed out in [7J. The g-Airy function A q satisfies g-difference equation: 

(1.6) A q (z) = A q (qz) - qzA q (q 2 z), A q (0) = 1 

with power series expansion 

q k2 (-z) k 



:i.7) aw = i + Ett 



k=l 



g )(l_g2)...(l_g 



fcV 



The purpose of this paper is to justify the above procedure for a more general class of 
orthogonal polynomials. More precisely we consider the family of orthogonal polynomials 
{p n (x, c)} defined by the recurrence relations 

(1.8) 2xp n (x, c) = p n+1 (x, c) + q~ nc (3 n (q, c)p n „i(x, c), 

where c > and f3 n (q,c) = 1 + o(l) as n — > oo. We will also remark on the case when 
(3 n (q,c) = 1 + YlJLi a>j(c)q X: > n for a sequence {A^} with A J+ i > Xj > 0,j = 1,2,3,..., 
in Section 5. This generalizes the g _1 -Hermite model. A second result generalizes the 
Stieltjes-Wigert and g-Laguerre models. We consider the three term recurrence relation 

(1.9) xq 2n+a+l p n (x) = a n p n+1 (x) + b n p n _i(x) + c n p n (x), 
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where, modulo linear scaling of x by ax + 5, 

oo oo oo 

(1.10) a n = l + Y,an, k q akn , b n = 1 + $^& n , fc g An , c n = 1 + ^ c n , fc g^\ 

fc=i fe=i fe=i 

where > 0, (3k > 0, and 7fc > for all k > 0. We will see that this recurrence includes 
both the g-Laguerre and Stieltjes-Wigert polynomials as special cases in next section 
where we will also introduce a one parameter family of solutions to the g-difference 
equation in (jl.fip . 



2 Main Results 

As in the case of f 1 1.11) we transform the polynomials p n (x,c) to the functions s n (t) 
defined as s n (t) = q n ^ 2 t~ n p n (x n (t),c) with x n (t) given by (jl.2p . Thus S-i(t) = 0, 
so(t) = 1, and 

(2.1) (l - Yj s n (t) = s n+l {q l lH) + i c (g, n)^^?- 1 / 2 *), n > 0, 

where \c(q, n)\ < K and, for each q, limn^oo c(q, n) = 1. 

Note that {s n (t)} are polynomials in 1/t 2 of degree / 2 J . 

Our main result is the proof of the convergence of {s n (t)}. We remark that the 
recurrence (12. ip indeed includes the recurrence relation of the g _1 -Hermite (II. 3p as a 
special case. 



Theorem 2.1. Let {s n (t)} be defined by the recurrence given by (12. ip with s_i(t) = 
and so(t) = 1. The limit lim n _ s>00 s n (t) exists and the convergence is locally uniform for 
teC\{0}. 

By the discussion immediately following (II. 3p . we have the following corollary. 
Corollary 2.2. Let {s n (t)} be defined by the recurrence given by (12. ip . We have 

lim s n (t) = Ag(^) 

locally uniform for t £ C \ {0}. 

The main ideas used in the proof of Theorem 12. II for the treatment of the family that 
generalizes the g _1 -Hermite polynomials can be adapted to the family of recurrence re- 
lations (i.e., ( II. 9p ) which include g-Laguerre and Stieltjes-Wigert polynomials as special 
cases. Due to the nature of non-symmetry, the implementation of these ideas become 
lengthier in the new situation. To transform (II. 9p . set 

x = Xn (t) = q- 2n ~ a t 
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and define S n (t) = q n (— t) n p n (x n (t)). Then the recurrence relation (II. 9p can be written 

as 

(2.2) S n {t) = (l-q n+l )S n+l (qh) + (l + q)c(q,n)(qt)- l S n (t) + qd(q^ 

where \c(q,n)\ < K and \d(q,n)\ < K 2 with lim^oo c(g, n) = lim^oo d(g, n) = 1. It is 
clear that S n (t) is a polynomial in 1/t with 

(2.3) S n (oo) = l/(g;g) n . 

To see that g-Laguerre polynomials {L^n\x; q)} and Stieltjes-Wigert polynomials 
{S n (x;q)} are special cases of (11. 9p and (ll.lOp . we recall that {L^\x;q)} is generated 
by 

(2 4) -xq 2n+a+1 L^(x; q) = (1 - g" +1 )4t(^ ?) + ?(1 - 9^)41 (*! ?) 

-[(i-g n+1 ) + g(i-g" +Q )]4 a) 

with initial conditions 

(&; ?) = !, m (&; = ; 

while the Stieltjes-Wigert polynomials {S n (x;q)} are defined by 

(2.5) - xq 2n+1 S n (x; q) = (1 - g" +1 )S n+1 (x; g) - [1 + g - q n+1 }S n (x; q) + g^ n _i(x; q), 
with the initial conditions 

S (x; q) = 1, Si(x; g) = ^ — — . 

1 - q 

Recall that (a; g)^ = (1 - a)(l - ag) • • • (1 - ag™" 1 ) • ■ ■ = U^oi 1 ~ af 1 )- 

Theorem 2.3. Let {S n (t)} be defined by (12.21) with initial conditions So(t) = 1 and 
Si(t) = [1 — g(l + g)c(g, l)]/t)/(l — g). Then the limit lim^oo S n (t) exists and 

—A (- 

{q;q)oo t 

locally uniform for t e C \ {0}. 

We have seen that y = A q (z) satisfies the functional equation 

(2.6) y{z) - y(qz) + qzy(q 2 z) = 0. 
It is worth noting that the family of functions 



lim S n {t) = ^A q {-), 



oc 



(2.7) F q (z; s) := £ 

„=-oo ^' q)n+s 

are solutions to (12.61) for all s. The notation used in the above equation is (cf. [6j IB]) 
(2-8) (a;g) fc = 7 %^. 

Observe that F q (z; s) is periodic in s of period 1, so there is no loss of generality in 
taking s in the strip Re(s) G [0, 1). Indeed, A q (z) = F q (z; 0). 



5 



3 Proof of Theorem 12.1 



Let {s n (t)} be given by the recurrence relation (12. ip with s_i(£) = and s (t) = 1. We 
will need the following lemma for the proof of Theorem 12.11 

Lemma 3.1. // there exist constants p > and M > such that 

(3.1) max\s n (t)\ < M 

\t\=p 

for all n > 0, then {s n (t)} is a normal family on C \ {0}. 

Proof. Since s n (t) are polynomials in 1/t 2 , they are all analytic on C\ {0}. So, by the 
maximal modulus principle for analytic functions, ( 13. ip implies that 

(3.2) max \s n (t)\ < M 

\t\>p 

for all n > 0. Next, we show that f]3 .21) implies that 

(3.3) max \s n (t)\ < (1 + ^^)M 

for all n > 2. 

Indeed, from the recurrence relation (12.11) . 

s n+1 (q 1/2 t) = ^1 - s n (t) - ic^Ti)^^!^- 1 / 2 *). 

Thus, using (13. 2p . for |t| = p, we have 

l^ + i(g 1/2 t)|< + M + W M= i 1+1 ~W~) M ' 

for n > 1. This proves (13. 3p . Repeating the argument above, we obtain that, for 
k = 1, 2, we have 

-I I TV- 

(3.4) max |s n (t)| < M TT(1 + ) 
for n > k + 1. 

Since limfc^oo g fc//2 p = 0, we see that the sequence {s n (t)} is uniformly bounded on 
compact subsets of C \ {0}. □ 

We need one more auxiliary result for the proof of Theorem 12.11 

Lemma 3.2. We have, for \t\ > 1, 

n K 
(3-5) Mt)\<l[(l + q k )A q (- w ) 

k=0 

for all n > 0. 
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Proof. We use induction. First it is trivial to verify that (13. 5p holds for n = 0. Now, 
assume (I3.5P is true for n and let us verify that (13. 5p is also true when n is replaced by 

From (12. ip . we have 



s n +i(t) = (1 - -p-)s n (q- 1/2 t) - -c(n,q)s n ^ l (q~H). 



Thus, for |t| > 1, 



< (n u n+i )Ki<r h2 t)\ + ^\s^i(q-H)\ 



\t\ 

< (i + na + ?*) Ai(-^i) + nV + ^c- S 

fc=0 ' ' ' ' fc=0 ' ' 

< 



n(i + /)U(-S + fiA(-S) 

i — n V I i II 11/ 



fc=0 
ra+1 



n( 1 + 9*)Ar(-|§. 

fc=0 



where in the last equality, we have applied (I1.6P with z = — K/\t\ 2 . This completes the 
proof of the proposition. □ 

Now, we are ready for the proof of Theorem 12.11 

Proof of Theorem \2.1[ First, note that (i) the infinite product in (13. 5p converges as n — > 
oo, (ii) A q {— M[) > 0, and (iii) when \t\ = 1, 



for some A > 0. Hence, with the help of Lemmas 13. II and 13. 2\ we see that {s n (t)} is a 
normal family for t G C \ {0}. 

To establish the convergence of {s n (t)}, we study the coefficients of s n (t). Recall 
that s n (t) is a polynomial of degree [n/2\ in 1/t 2 . Write 

( 3 -6) S n (t) = ^2 a n,k^- 

k=0 

Then, the fact that {s n (t)} is a normal family implies that {a n) k\ is a bounded set. 
Assume that 

(3.7) \a n , k \<M, for < k < [n/2\, n = 0,1,2,.... 
Next, we show that, for every fixed k, 

(3.8) lim a Ujk =: f k 



exists. 

For n > 0, using (13. 6p in ( 12. ip . we get 

k=0 fc=0 y fc=0 

Comparing the constant terms on both sides of the above equation, we obtain 

a n,0 — a n+l,0i n — 0, 1, 2, .... 

From s (t) = 1, we get a ,o = 1, so, we have 

(3.9) a„,o = l, 71 = 0,1,2,.... 

Next, for k > 0, comparing the coefficients of l/t 2k , we get, for large n, 

n a n+l,k k-l 1 \ 

a n ,k - q a n ,k-\ = z h ? c(g, n)a n _i jfc _i 

Or, equivalently, for large n, 



(3.10) q k a n>k - q n+k a n ^_ x = a n+hk + q 2k 1 c(q, n)a n _ iife _ 



Replacing n by n + I in ( 13.101) , we get 

(3.11) q k a n+ i )k - q n+l+k a n+ i tk -i = a n+i+ltk + g 2fc_1 c(g, n + Z) a n+«-i,fc-i- 
Subtracting ( 13.111) from ( 13.101) and re-arranging the terms, we get, for large n, 

k ( \ tl~ \~ k ( I \ 

0>n+l,k — On+J+1,* = q \ a n,k ~ Q>n+l,k) ~ Q \P"n,k-l ~ Q a n+l,k~l) 

- q 2k ~ l [c{q, n)a n _ life _! - c(q, n + Z)a n+J _i )fc _i]. 
This, together with (13. 7p and the fact that < q < 1, gives us 



\ a n+l,k ~ a n+l+l,k\ 



< q k \ a n,k — 0"n+l,k\ + q n+k 2M + q 2k 1 |a n _ifc_l — CL n +l-l,k-l \ 



(3.12) +q 2k ~ 1 M(\c(q,n) - 1| + \c(q,n + I) - 1|). 

With the help of (13.121) . we are ready to prove (13.81) by induction on k. 
When k — 0, equations in (I3.9P give us 

(3.13) lim a n fi = 1. 

Now, assume ( I3.8P is true when k is replaced by k — 1. Then, for any £ > 0, there 
exists a number iV > such that 

(3.14) |a n -i,fc-i - a n+i _i ifc _i| < e, n> N, I > 0. 
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From lirrin^oo c(q, n) = 1, we may choose N large enough to ensure that 

(3.15) \c(q,n) - 1| < e, n > N. 
Thus, for n > N, (13TT2D implies 

(3.16) |a n+1 , fc - a n+l+1>k \ < q k \a n , k - a n+Lk \ + q n+k 2M + q^e^l + 2M). 
After repeatedly using (13.1 6p m times, we arrive at 

\n n I <~ n( m + 1 ) k \n n I 

\ a n+l,k ~ a n+l+l,k\ — y \ a n-m,k ~ a n-m+l,k\ 

m m 

+ q n+k 2M (l j{k ~ l) + q 2k ' l e{^ + 2M) J2 V 3k > 

3=0 3=0 

< q {m+l)k 2M + q n+k 2M ^—r + q 2k ~ l e(l + 2M)—^, n-m>N, I > 0. 

Take m = [y/n\. Then lim^oom = oo and lim ri ^ 00 (n — m) = oo. Therefore, for large 
n and for all I > 0, 

q2k-i 

lim \a nk - a n+l)k \ < e(l + 2M)- r. 

n— >oo 1 — q K 

So, 

(3.17) lim (on.fe - a n+itk ) = 

n— >oo 

uniformly in I > 0. Thus, {a n ,fc}n>fc is a Cauchy sequence. Hence, (I3.8P is true for any 
fe = 0,1,2, ... 

Finally, we need to derive the convergence of {s n (z)}. Although this can be done 
by a bounded convergence argument based on ( 13. 7ft and ( 13. 8ft . we choose the following 
more elementary argument for its simplicity: Note that every subsequence {s n (t)} n( =A of 
{s n (t)} has a sub-subsequence {s n (t)} neAl (A 1 C A) that converges locally uniformly in 
C\ {0}. Assume 

f{t):= lim s n {t). 

Then 

/W(f) = lim s g>(t). 

n— >oo,nGAi 

Evaluating at t = oo, we see that 



oo ^ 
/(*) = 5^/*72fc- 



3=0 

Therefore, the whole sequence {s n (t)} must converge to the same function /(i) locally 
uniformly in C \ {0}. □ 
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4 Proof of Theorem 12.3 

We need to recall some notations from g-theory and solve a g-difference equation. 



4.1 The Functional Equation 

Recall that (see <I2T5|> ) 

n-l 

(4.1) (o;g) = l, (a; q) n = JJ(1 - ag J ), 

3=0 

and we note the Euler identity P, Page xvi, (18)] 

00 2 n 1 

We will need to solve the functional equation 

(4.3) f {z ) = f( zq *) + ^f {zq -*) + a l±A f ( z)] a = ± i. 

z A qz 

Proposition 4.1. Up to a scaling factor, the above functional equation has a solution: 

(i) when a = 1, f(z) = A q (-), and 

(ii) when a = —1, f(l/z) is an entire function with positive Taylor series coefficients. 

Proof. Assume that f(z) = Yl^Lo fn z ~ n ■ Since /o 7^ is a scaling factor there is no loss 
of generality in assuming / = 1. By equating coefficients of z~ n on both sides of (14. 3 p 
we see that 

(4.4) a(l-r 2n )/n = 9 2n -7n-2 + (l + ?)?- 1 /n-i, n > 0, 

with fi = aq/(q — 1). Now, in order to normalize the coefficients in difference equation 
(@~3D above, set f n = g n q n2 . Thus g = 1, g x = a/ (q - 1) and 

g n (q 2n - 1) = aqg n - 2 + a(l + q)g n -i, n>0. 

We now solve the above recursion using generating functions, so we set G(z) = Yl^Lo 9nZ n - 
In view of the initial conditions, the functional equation for G(z) is 

(4.5) G(z) = G(q 2 z)/[1 + a(l + q)z + aqz 2 ]. 

If a = 1 then the denominator in (14. 5 p factors as (1 + z)(l + qz), and iterating (14. 5 p with 
the initial condition G(0) = 1 gives 
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by the Euler identity. In the last step we used (I4.2p . This leads to f(z) = A q (l/z). 
On the other hand when a = — 1, the denominator in (14. 5 p factors as 

1 - (1 + q)z - qz 2 = (1 - z/a)(l - z//3), with a < < (5 and /3 < \a\. 

In this case 

G(z) ' 



z/a; g 2 )oo(2//3; 



Since < (3 < \a\, we see that the singularity of G with the smallest modulus is 0. 
Applying Darboux's asymptotic method, [12], we conclude that 



(3~ n 

3 " = (,WUP/a;^ [1 + 0(1)l ' 



/ OO 



Thus /„ = (7nQ ,n and /(I/2;) is an entire function of z. Moreover the recurrence relation 
( 14. 4p and the initial conditions show that f n > 0. Therefore the function f(l/z) is an 
entire function of z with positive Taylor series coefficients. □ 

Proof of Theorem \2.3l We follow the ideas used in the proof of Theorem 12.11 We will 
skip similar argument. The main task is to show that the family {S n (t)} defined by 
( 12.21) is a normal family for t G C \ {0}. 

Let f b {t) be the solution with f b (oo) = 1 in Proposition 14. II when a = —1. We claim 
that for \t\ > 1, 

(4.6) m) \< ™* 



nLi(i-^)' 

for n = 0,1,2,... 

Let us use induction to verify (14. 6p for n = 0, 1, 2, ... Clearly, from the initial condi- 
tions and the fact that f b {\t\/K) = 1 + Kq/((1 — g)|t|) + (positive terms), we see that 
(14.61) is true when n — 0,1. Now assume that (14.61) is true for n. 

From (I2.2p . we have 

(1 - q n+l )S n+1 (t) = S n (q- 2 t) - (1 + q)c{q,n){q-H)- l S n {q~H) - q 5 d(q,n)r 2 S n ^{q-H). 
So, for \t\ > 1, 

\(l-q n+1 )S n+1 (t)\ < \ Sn (q-H)\ + { -^^\S n (q-H)\ + ^\S n ^(q-H)\ 



t 



< f b (q- 2 \t\/K) f | K(l + q)q \ | K\* f\q+\t\/K) 



< 







-q k ) 




1 




nLi 


(i 


-q k ) 




1 






(i 


-q>) 



x + K(l + q)q \ /6( g- 2 jt| + K^ fb/ q-*\t\ 



t J K t 2 J v K 



f\- 
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where in the last step, we used the functional equation ( f)4.3p with a = — 1) satisfied by 
f b with z = q~ 2 \t\/K and the fact that f b {\) has positive Taylor series coefficients. This 
implies (14. 6[) holds when n is replaced by n + 1. Therefore, by induction, (14.61) is true 
for all n = 0, 1, 2, ... 

The rest of the proof goes like the one for Theorem 12.11 From the claim, we can 
easily show that {S n (t)} is a normal family for t G C \ {0} and from here and working 
with the recurrence relations, we can show (by a similar argument as given in the proof 
of Theorem 12 .ip that the coefficients of S n (t) converge: if 

n 

Sn(t) = a n,kt k 

k=0 

then 

lim a n)k = f k , for some f k , k = 0, 1, 2, 

n— >oo 

which immediately yields the convergence of the whole sequence {S n (t)}. 

Finally, writing the limit as f(t): f(t) : = lim^oo S n (t), and let n — > oo on both sides 
of (ED to get 

/(*) = ftft) + (1 + qXqty'fit) + qt- 2 f(q- 2 t). 
From this, with the fact that lim^oo S n (oo) = l/(g;g)oo, Proposition 14.11 when a = 1 
implies that f(t) = A q (\)/{q] q)^. □ 

5 Error Terms 

Unlike the asymptotics of the g _1 -Hermite polynomials [7], the terms beyond the main 
term in the asymptotic expansion of s n (t) seem to be intricate. 
In (12. II) we assume that c(q, n) has the convergent expansion 



(5.1) c(q,n) = Y,c k q nk + o(q nm ), c = 1. 

fc=0 

We further assume that s n (t) has the asymptotic expansion 

m 

(5.2) s n (t) = fnA±/t 2 )q nk + o(q nm ). 



k=0 



Now substitute for s n from (15.21) in (12.11) and equate the coefficients of q nk . When k = 
we see that f n ^(l/t 2 ) solves 

(5.3) y(u)-y(u/q)-uy(qu) = 0, 

with u = 1/t 2 . When k > we conclude that 

fn,k( u ) ~ q k fn+i,k(u/q) - uq~ k f n -i,k(qu) 
(5 4) k ~ 1 

= w/n,fc-l(w) + U 2j Ck-jq~ 3 fn-\,Mu), 
j=0 
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for k = 1,2, ■ ■ ■ . Since s n (t) is a polynomial in 1/t 2 we expect f n ,k( u ) to be analytic in 
u in a neighborhood of u = 0. Thus f n ,o(u) = A q {u). Note that ( 15. 4 p implies / n ,fe(0) = 
for fc > 0, since s n (oo) = 1 = A 9 (0). Let g n ,k(u) = f n ,k(u)/u, for k > 0. Thus (15. 4p gives 

(5.5) SVi,i(u) - g n +i,\{u/q) - ug n - hl (qu) = A q (u) + cxA q (qu). 

The case c(g, n) = 1 — q n is very exceptional. In this case Ismail j7] showed that 

In this case ( 15. 4 p becomes 

, r _x fn,k( u ) ~ Q k fn+l,k(u/q) - uq~ k fn-Lkiqu) 

(5-7J , , 

= Ufn,k-l(u) ~ Uq fn~l,k-l{qu), 

and by induction f n ,k{u) must have the form u k g n ^{u) and g nj fc satisfy 

(5.8) #n,i(-u) - g n +i,\{ujq) - ug n - X)1 {qu) = g n ,k-i{u) - g n -i,k-i(qu), 

which leads to the solution given by (15. 6p . 

The analysis of (15. 5p in general seems to be complicated. As an illustration we first 
consider the case c(q,n) = 1 + c\q n . We then let g n> x = YljLoCn,j u ^ an d substitute it 
in ( 15.5ft . Thus Cn+i,o — Cn,o + 1 + ci, and we conclude that C„,o = n (o,o + n(l + c\). 
Therefore Cn,o is independent of n if and only if c\ = — 1. When ci = -1 we are led to 
the expansion (15. 6p . 

For general ci we make the Ansatz g n ,i(u) = nF{u) + G(tt) + o(l) which leads to 
F(m) = A q (u) and that G(w) solves 

(5.9) G(u) - G(u/q) - uG(qu) = 2A q (u/q) + c^qu). 

Writing G(u) = YIJLq ^j u ^ we nn d the following two-term recurrence relation for A n 

(5.10) ( i yxjq' 1 = i w V ; 1)2 a, . + 2 (-iyi±^, 

which is easy to solve. This process can be iterated but it gets complicated as we proceed 
to higher k. 

The error analysis for (12.21) is even harder. Ismail [7] proved 

(5.11) S n (q- 2n t; q) = — — £ Al^gtt V^( g -/t), 



and 



(5.12) 



<rL^{x n {t)-q) 1 ^ g^ 2 ,,,, 



x 



(-*)" (g; g)oo ^ (g; g) 7 

E {q'->q)m i ~i \m— s n sa+(m— 2s) 2 / 2 A I 2s-m 
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Since S n (oo) = l/(g;g)„ we let S n (t) = Y,T=o fn,k{ u )q nk wnere u = V*- Equa- 
tion (E2D leads to 

(5.13) f n ,o( u ) = fn+i.o{u/q 2 ) + 1 + Q uf nfi (u) + qu 2 f n -i fl (q 2 u). 

Q 

Any constant times A q (u) will solve (I5.13P as per (i) of Proposition 4.1. In general, we 
have 

fn,k(u) = fn+l,k(u/q 2 )q k + - ^Cfc-j/njM 

+gw 2 d k _jq~ 3 f n -i,j(uq 2 ) ~ u2 ^2 d k-i- 3 q~ 3 fn-l,j{uq 2 )- 

3=0 3=0 

We do not expect to depend on n. One reason is the following. In (I2.2p . let 
S n(t) = + Y!j=i an,jt~ J . We then have 



Thus 



a n ,i = (1 - <? n+ )on+i,ig + —, z-c(q,n). 

q{q; q)n 



2tj / \ 2r? 2/ \ — ^~ Q / \ In 

q {q;q)na n ,i = q (q; ?)n+i^n+i,i H c{q,n)q 

2n (n- n , = M 4-1 V"" 1 rfn ^«- 2 J 



which gives g 2 "(g; g) n a n ,i = (1 + 1/g) E?=o c (<?, j)<? 
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